One-Dimensional Extended States in Partially Disordered Planar Systems 
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We obtain analytically a continuum of one-dimensional ballistic extended states in a two- 
dimensional disordered system, which consists of compactly coupled random and pure square lattices. 
The extended states give a marginal metallic phase with finite conductivity <to = 2e 2 jh in a wide 
energy range, whose boundaries define the mobility edges of a first-order metal-insulator transition. 
We show current-voltage duality, Hu/T scaling of the conductivity in parallel magnetic field Hu and 
non-Fermi liquid properties when long-range electron-electron interactions are included. 
PACS numbers: 72.15. Rn, 71.30.+h, 73.20.Fz 



Anderson localization and the associated metal- 
insulator transition (MIT) remain one of the most ac- 
tive fields in condensed matter physics Q], The one- 
parameter scaling theory predicts localization of all the 
states in one or two dimensions for any amount of disor- 
der |], although even in one-dimension (ID) isolated ex- 
tended states can exist as a consequence of specific short- 
range disorder correlations HQ]. In a two-dimensional 
(2D) disordered system derealization is possible for non- 
interacting electrons only by breaking the time-reversal 
symmetry via a magnetic field p{ or the spin-rotation 
invariance by spin-orbit coupling ||. Extended states 
can also appear in continuous 2D models with random 
S- function potentials of infinitesimal action distance . 
The presence of extended states in 2D random systems is 
of great interest since metallic phase and MIT has been 
recently reported in Si-MOS planar structures p|-pO[ and 
GaAs-based materials pT| , in the absence of magnetic 
fields. 

We analytically show extended wave functions in par- 
tially random structure composed of two compactly cou- 
pled square lattices, one random and the other periodic. 
At the same time, according to the scaling theory the 
rest of the states are localized. Moreover, the extended 
states have special momentum and form a ID continuum 
which gives 2e 2 jh conductivity, in a finite energy range. 
The boundaries of this range define the mobility edges 
(ME) of a first-order MIT. These states are perfect (bal- 
listic) and coexist with localized states in their energy 
range. They correspond to a non-Fermi liquid, in agree- 
ment with recent scaling expectations for the presence 
of a MIT in 2D fjjj. We show that the system under 
long-range Coulomb interactions reduces to the marginal 
Fermi liquid proposed in Jll|. 

We consider two coupled square lattices with the sites 
of one located above the plaquette centers of the other. 
One plane is pure and the other is random. The tight- 



binding Hamiltonian reads 

H = ^ e\ 4 a\ z ax,i + ^ >\ : "l,,<i\., + H.c.) 

+ Yl *( a I,i°2j -+H.c.), (1) 

where is the destruction operator for an electron at 
site i of the A = 1(2) pure (random) plane, t\ is the 
nearest-neighbor (NN) hoppings within the A-th plane 
and the third term connects each random site to its 
NN pure sites via the hopping parameter t. We choose 
ti = t2 = 1, setting the energy unit throughout, although 
our conclusions are also valid if ti is random. The sites 
of the pure plane have energies e\j — and of the ran- 
dom plane e%,i, randomly distributed between —W/2 and 
W/2. For convenience, the length units and the origin of 
the coordinates are set so that the pure (random) sites 
count odd (even) numbers. 

Obviously, extended wave functions can have zero am- 
plitude on the random sites and finite amplitude on the 
pure sites, being not influenced by randomness. We can 
seek for such states in the considered partially disordered 
system. Indeed, for an L x L lattice with the described 
geometry and periodic boundary conditions in all direc- 
tions, with L even and x n , y n the coordinates of the 
n = (A, i)-th site, the states 



^x,k x { x n,y n ) = ^-sin(^Y-)exp(ik x x r . 



lp2,k y (Xn,y n ) 



;| sin (~Y") exp(ifc y y„), 



(2) 



(3) 



where k x (k v ) = 2j x (j y )ir/L, j x (j y ) = 1,2, ...,L/2, are ex- 
act eigenstates of the Hamiltonian H having zero am- 
plitude on the random sites. This can be easily verified 
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by applying H on ^i,k x ip2,k y (x n ,y n ) to obtain 
their eigenenergies 

E = 2[cos(2fc x ) - 1] and E = 2[cos(2fc y ) - 1], (4) 

respectively. The transverse momentum of these states 
is fixed to k y {k x ) — n/2 and their longitudinal momen- 
tum k x (k y ) runs in a II? Brillouin zone, forming a ID 
continuum in the energy range [—4, 0] [E3. 
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recently suggested, from experimental data in 2D hole 



systems of GaAs heterostructures |15 



The obtained conductivity a = <to is invariant under 
scaling so that from the viewpoint of the scaling theory 
it does not imply a true metallic phase for the system, 
but a critical phase instead, with the corresponding /3- 
function only reaching zero Moreover, since the ex- 
tended states coexist with localized states in a wide en- 
ergy range, having lower dimensionality and measure, we 
can rename the obtained phase "marginal metallic" , fol- 
lowing the terminology "marginal Fermi liquid" proposed 
by Varma el al (jlj. 
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FIG. 1. The density of states p(E), in arbitrary units, for 
a finite 2D random sublattice system of linear size L — 52 
with t — 1.2 and W = 1.5 by taking averages over 100 ran- 
dom configurations. The scale invariant conductivity ao is 
also shown. 

The perfectly extended states of Eqs. (||) and (||) 
have momentum along both the 2D principal axes and 
do not violate the scaling theory, since they effectively 
decouple from the random sites although the two sublat- 
tices are compactly coupled. At a given energy E, with 
—4 < E < 0, they provide propagating channels with 
(To = 2e 2 /h conductivity (if spin is included) along the x 
(or y) direction, independently of the system size L. In 
Fig. 1 we plot an and the density of states p(E) as a func- 
tion of E. The rest of the 2D states are asymptotically 
localized, as expected from the scaling theory [Q. The 
transparent channels give minimum metallic conductiv- 
ity ao = 2e 2 jh and one obtains a first-order MIT from 
an to at E c \ = —4 and E C 2 = 0, which can be re- 
garded as mobility edges. For the chosen structure we 
observe that the lower ME lies near the bottom of the 
band (see Fig. (1)), so that the Fermi energy can move 
into the metallic regime even for small electronic doping. 
For example, this is the case in silicon MOSFETs where 
the metallic phase is reached at very low electron densi- 
ties (~ 10 n cm- 2 ) |§-|l§. It should be also noted that 
a minimum metallic conductivity close to e 2 /h has been 
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FIG. 2. The current-voltage curves at several 

Fermi levels around the mobility edge with temperature 
k B T = 1.1 x 10~ 5 ti. 

The step-like conductivity a{E) = ^-0(E - E cX ) at 
E c \ gives nonlinearity of the current- voltage (I — V) 
curves with I — V duality. If the sample is connected to 
two metallic leads with the same Fermi level Ep, an ap- 
plied voltage V between the two contacts opens a trans- 
mission window of width eV around Ef at zero tem- 
perature, due to the electronic occupation in the leads. 
Nonlinearity occurs if the window boundary passes across 
the ME when changing the voltage. It becomes smoother 
at finite temperature T. In order to illustrate this effect 
we compute the current from the formula 



I(V) = 



In 



k B T 



E C2 



1 



1 + exp[(£ -E F - eV/2)/k B T] 
1 



1 + exp[(£ - E F + eV/2)/k B T] 



dE, 



(5) 



where In is a prefactor. In Fig. 2 the calculated I — V 
curves display nonlinearity if Ep ^ E c ±. For Ep larger 
(smaller) than E cl the resistivity increases (decreases) 
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with |V|, giving the symmetric behavior about Ep = E c \ 
known as I — V duality, exchanging the insulating with 
the metallic phases. Similar kind of duality has been 
observed in silicon MOSFETs (Fig. 1 of Ref. JuJ). 

The extended electrons, decoupled from the rest, form 
a pseudo-Fermi sea embedded in the localized electrons. 
The pseudo-Fermi surface has only four points corre- 
sponding to two Luttinger liquids crossing with each 
other. This provides a possibility for applying the ID 
or quasi-l-D Luttinger liquid theory to 2D systems, as 
pointed out several years ago p6| |. Nevertheless, the sys- 
tem should be regarded as 2D non-Fermi liquid owing to 
its measure, its ID-like spectrum, as well as the presence 
of the background from the localized electrons. 

If Coulomb interactions are included the corresponding 
matrix elements read 

(/ii,/i 2 |F|M3,M4> = c i( n i) c 2( n 2)^(|ni - n 2 |) 

ni,n 2 

c 3 (ni)c 4 (n 2 ), (6) 

where /Zj is a single-particle state of H and Cj(ny) its coef- 
ficient at site r\j . We are interested in elements where one 
of the initial states (/is) is extended and the other (/X4) is 
localized which describe scattering between extended and 
localized electrons. The extended states have amplitudes 
with rapidly changing sign in one direction, while the lo- 
calized states show no phase coherence having randomly 
varying amplitudes. Thus, for long-range Coulomb inter- 
actions the elements with localized final states /ii and /i 2 
are almost zero, due to cancellations arising from the pos- 
itive and negative signs of the extended waves. This leads 
to an extremely small probability for an extended state 
to be scattered into a localized state, although the num- 
ber of localized states is much higher Jl4| . If a final state 
(e.g., /^i) is also extended with the same wave symmetry 
as the initial state the matrix element is of order ~ 1, 
since the product cj(ni)c 3 (ni) eliminates the sign of the 
extended wave coefficients. The above discussion implies 
that the extended states found can survive in the pres- 
ence of long-range Coulomb interactions. However, they 
experience fluctuations from particle-hole excitations of 
the localized electrons, which can be described by the 
polarizability 



2k B T 



E 



G (ei, iv n )G (e 3 , iv n + iu m ), (7) 



where Go is the single-particle Green function of the lo- 
calized states, the sums for i and j are over the relevant 
localized states for the given interaction, ti is the energy 
of state i, and v n = (2n + VjitksT with n an integer. 
After summing over the Matsubara frequencies and as- 
suming constant density of the localized states, one has 



Imll(cj) 



A {w/k B T), for |w| < k B T, 
Apsignoj, for \oj\ > k B T 



where Aq is a prefactor dependent of density of states. 
This is the basic hypothesis of the marginal Fermi liq- 
uid theory |l3). The key point in this derivation is the 
momentum independence of the Green function for the 
localized states. If there are no added impurities in the 
host sublattice, this manifold is a pure marginal Fermi 
liquid which can give finite conductivity at zero temper- 
ature, distinct from the non-pure marginal Fermi liquid 
which experiences further impurity scatterings |17[j. 
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FIG. 3. Main figure: the H\\/T scaling of the conduc- 
tivity change due to a parallel magnetic field for the system 
at the transition point. Inset: the Fermi energy Ef vs. (f>, 
the flux through a primary inter-plane triangle with the field 
along the x axis, averaged over 100 random samples with size 
L = 30 (2 x 15 x 15), and the rest of parameters as in Fig. 1. 
The units shown are derived for a sample with a = b — 6 A, 
t\ — leV. The slope of the inset curve gives gi — —\.\\ib- 



An in-plane magnetic field applied to the quasi-2Z3 
system in the Coulomb gauge gives in-plane NN hop- 



t 



(8) 



exp(i</> shitty) and = exp(— i^sinccy), 
where cty is the angle between the hopping bond ij 
and the field, <f> — hH\\ab is the flux through a pri- 
mary inter-plane triangle, a the planar lattice spacing 
and b the distance between the planes. The inter-plane 
NN hoppings are unchanged in this gauge. Thus, the 
states of Eqs. (||) and (||) survive with eigenenergies 
E = 2[cos(2fc2; + 2(/)sina a ;) — cos(2cf>s'may)) and E = 
2[cos(2/c 1 ,-|-2(/)sinQ;j ; )— cos(2(/)sin a x )], respectively, where 
a x ( y j is the angle between the x(y) axis and the field. It 
can be seen that the shift of E cl (or E c2 ) is of the or- 
der of H| for small fields. At the same time, the whole 
band including the localized states changes, so that the 
Fermi level Ep is also shifted if the electron density is 
kept fixed. In order to obtain the field dependence of 
Ep we compute the energy spectrum of a finite random 
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system and the result is shown in the inset of Fig. 3. In 
general, A(Hu), the energy difference between E c i(H\\) 
and Ep (H\\), can be expanded as a power series of the 
field, via A(iT||) = A(0) + gxH\\ + .... For small fields 
we can only keep the linear term. On the other hand, 
the spin degeneracy is removed by ±[1bH\\ for up and 
down spins, with fis being the Bohr magneton. For a 
small bias voltage the conductivity can be expressed in 
terms of Landauer-Biitticker formula Q and at the MIT 
(A(0) = 0) one has 



a(H h T) - a(0,T) 



E 



1 



h — r i 

g=±l t 



exp 



(gi+gM-B)g|| 

k B T 



(9) 



which shows H\\jT scaling behavior. In Fig. 3 the H^/T 
dependence of conductivity is plotted for one sample, 
showing good agreement with Fig. 3 of Ref. jl9| ) . 

In summary, we have rigorously shown a band of ID 
ballistic extended states which propagate freely in a 2D 
mixed pure-random system, without magnetic field or 
spin-orbit coupling. The ID extended continuum is lo- 
cated at the edges of the first Brillouin zone (k x = ir/2 
and any k y or k y = n/2 and any k x ) of the 2D square 
lattice. Moreover, the ID extended states coexist with 
localized states, in a range of energies whose boundaries 
define a first-order MIT. We derive I — V duality, H\\ /T 
scaling, and demonstrate non-Fermi or marginal Fermi 
liquid behavior. These features have not been obtained 
previously in the presence of disorder and are of pro- 
found significance in understanding many experimental 
findings, especially the 2D MIT and the underdoped 
cuprates. In this Letter we only focus on the general 
characteristics of a class of partially disordered materi- 
als, with at least one sublattice perfect. In Si MOSFETs 
and GaAs-AlGaAs heterostructures, electrons or holes 
are confined within several atomic planes and due to lat- 
tice mismatch, or atomic diffusion, the planes adjacent to 
the interface are random, while some planes relatively far 
from the interface are less random and may support bal- 
listic channels. Another example presents YBa2Cu306+a; 
where the Cu02 planes are perfect and the CuO^ planes 
random due to oxygen vacancies. Although for these ma- 
terials more realistic models are clearly needed, possibly 
by including more layers, our results may still shed light 
on many of their common features. 
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